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ABSTRACT. General Relativity gives that finitely many point masses between an observer and a 
light source create many images of the light source. Positions of these images are solutions of 
r(z) = z, where r(z) is a rational function. We study the number of solutions to p(z) = z and 
r(z) = z, where p(z) and r(z) are polynomials and rational functions, respectively. Upper and 
lower bounds were previously obtained by Khavinson-Swiatek, Khavinson-Neumann, and Petters. 
Between these bounds, we show that any number of simple zeros allowed by the Argument Principle 
occurs and nothing else occurs, off of a proper real algebraic set. If r(z) = z describes an n-point 
gravitational lens, we determine the possible numbers of generic images. 



1. Introduction 

One of the results of Einstein's General Theory of Relativity is that a point mass placed between 
an observer and a light source will create two images of the source. If this single mass is replaced 
with a distribution of masses, significantly more complicated configurations of images can be cre- 
ated. Multiple images were first observed by astronomers in the 1970's and further technological 
advancements pushed gravitational lensing as an important tool in astrophysics. Gravitational lens- 
ing has also become an exciting field of research in mathematical physics — see the recent beautiful 
surveys [|5] [8] [12l [T4J and for deeper discussion, including the history, see [fT3l and |[T6l . 

Suppose that the distribution of mass is well-localized relative to the distances between the 
observer and the mass and relative to the distances between the masses and the light source. Images 
of the light source are described by solutions for z to 

(...) 

where fi is a compactly supported measure describing the distribution of mass projected onto the 
plane through the center of mass perpendicular to the line from the observer to the light source. 
See, for example, [21 and tfTTl . 

An early and important result in gravitational lensing, due to Burke ||2), is that if fx is a smooth 
mass distribution, then the number of solutions to (ll.lt is odd. This was generalized by Petters 
[fTT| to the situation where \i has smooth density except at g points. In this case, he showed that 
the number of solutions is congruent to (g — 1) modulo 2. See also Ifl4l Thm. 1]. 
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In this paper, we will focus on the case of n point masses. For 1 < j < n, let Uj be a positive 

n 

mass located at Zj. In this case, \i = aj5(z — Zj) and (11.11 ) simplifies to become 




(1.2) 



For the remainder of the paper, we refer to (|1.2t as the lens equation. 

It is also interesting to study variations of Equation (11.21) . replacing the sum on the right-hand 
side with a general polynomial p(z) or rational functions r(z). In the polynomial case, Khavin- 
son and Swiatek [9] used a clever application of the Fatou Lemma from holomorphic dynamics 
combined with the "Argument Principle" [31 to show that if p(z) has degree n, then the number 
of solutions to p(z) = z is bounded above by 3n — 2. It was a delicate question of whether this 
upper bound was achieved for each n. Using Thurston's Theorem from rational dynamics, Geyer 
[6] proved the sharpness of this bound. 

In the rational case, Khavinson and Neumann Q used similar techniques as to prove that 
if r(z) has degree n, then the number of solutions is bounded above by 5n — 5. (Note that 
r(z) = p(z)/q(z) has degree n = max(degp, degg).) Surprisingly, the sharpness of this bound 
for each n had already been proved by Rhie 03], using an explicit construction of an appropriate 
configuration of masses in the lens equation (I1.2I ). 

Let U be an open subset of C. A function / : U — > C is called harmonic if both the real and 
imaginary parts of / are harmonic in the classical sense. A zero z = x + iy of f(x,y) = 
u(x, y) + iv(x, y) is called simple if the Jacobian T>(z ) = det | | ^ 0. We call a polynomial 
p(z) simple if all of the zeros of p(z) — z are simple. Similarly, we call a rational function r(z) 
simple if all of the zeros of r(z) — z are simple. 

If p(z) is a simple polynomial of degp > 2, then there is also an obvious lower bound on the 
number of solutions by degp, as a consequence of the "Argument Principle." If r{z) = p(z)/q(z) 
is a simple rational function of deg r > 2, then there is a lower bound on the number of solutions 
depending on the degrees of p and q. For example, if degp < degg, then the lower bound is 
deg r — 1 . 

A special subcase of the rational case is obtained by considering rational functions of the form 
(11.21) with all positive masses a-;. Since Rhie's examples were constructed with positive masses, 
the upper bound of 5n — 5 is still achieved. Meanwhile, Petters [11] showed using Morse Theory 
that if r(z) is a simple rational function of the form (11.21 ), then the lower bound on the number of 
images is n + 1. 

We will look at each of the three cases mentioned above (polynomial, rational, and physical) 
from the perspective of "parameter spaces" and with a motivation of understanding what numbers 
of solutions between the lower and upper bounds can occur generically. 

For the purposes of this paper, we will parameterize the space of polynomials by their coeffi- 
cients, letting 

V n = {a n z n + a n _ x z n ~ x + ■ ■ ■ + a | a n E C \ {0} , aj E C for < j < n - 1} ^ C \ {0} x C n . 



Remark. After performing a rotation and/or shift in z, one can suppose that a n E IR + and a n _i = 
without affecting any of the statements below. However, we consider the present definition of V n 
more natural. 
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FIGURE 1 . Parameter space pictures for slices of A4 2 , M3, and M4, respectively. 
The fixed masses are denoted with large, black dots. 

Left: There is one fixed mass located at (1, 0). Placing a second, equal mass in the 
green region will produce 3 images and placing in the red region will produce 5 
images. 

Center: There are two equal, fixed masses located at ^— |, and ^— |, — ^ 

Placing a third, equal mass in the green region will produce 4 images, red- 6 images, 
blue- 8 images, and yellow- 10 images. 

Right: There are three equal, fixed masses located at ^— |, , ^— |, — , and 

(1,0). Placing a fourth, equal mass in the magenta region will produce 1 1 images, 
light blue- 13 images, and black- 15 images. 



Theorem 1.1. Let SV n (k) be the set of simple polynomials of degree n > 2 with p(z) — z having 
k roots. Then SV n (k) is a non-empty, open subset ofV n if and only ifk = n, n + 2, . . . , 3n — 2. 
Furthermore, the set of non-simple polynomials, NV n , is the complement of the union 

U SV n (k) 

fc=n,n+2,...,3n— 2 

and is contained in a proper real algebraic (hence measure 0) subset ofV n . 

We parameterize the space of rational functions of degree n by their coefficients up to scalings, 
with the condition that p and q are relatively prime. More specifically, let 

a n z n + . . . + a 



r(z) 



(1.3) 



b n z n + . . . + 60 

for (a n , . . . , do, b n , . . . , b ) G C 2n+2 . The space 1Z n of rational functions of degree n can be param- 
eterized by 2n + 2-tuples of complex numbers (a n , . . . , a , b n , . . . , b ), considered up to non-zero 
complex scaling with two restrictions: 

(1) a n 7^ or b n ^ and 

(2) The resultant of a n z n + . . . + a and b n z n + . . . + b is not equal to 0. 
As such, it is an open subset of CP 2n+1 . 

Theorem 1.2. Let STZ n (k) be the set of simple rational functions r(z) of degree n > 2 such that 
f(z) = r(z) — z has k roots. Then for all n > 2, S1Z n (k) is a non-empty, open subset oflZ n if and 
only ifk — n — l,n + l,...,5n — 5. Moreover, the complement of 
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U sn n (k), (1.4) 

k=n— l,n+l,...,5n— 5 

is contained in a proper real algebraic (hence measure 0) subset of lZ n . 

Remark. In the rational case, the complement of the union given by Equation (11.41 ) consists of both 
the hyperplane b n = and the set of non-simple rational functions, NlZ n . 

In the physical case, we parameterize the space of all configurations of n positive masses in C 
by n-tuples 

M n = {{{zi, <7i), . . . , {z n , a n )) e (C x M+) n | Zi ^ z d if i ^ j}. 

Note that our parameterization represents the masses as "marked," i.e. if two masses have the 
same mass and are interchanged, then the corresponding point in Ai n is different even though the 
physical configuration is the same. 

Theorem 1.3. Let SM. n {k) be the set of simple, "positive massed" rational functions r(z) G M. n 
having degree n > 2 that yield k solutions to Equation (11.21) . Then for all n > 2, SM. n (k) is a 
non-empty, open subset of \M n if and only ifk = n + 1, n + 3, . . . , 5n — 5. Moreover, the set of 
non-simple positive massed rational functions, NA4 n , is the complement of the union 

U SM n {k) 

fc=n+l,n+3,...,5n— 5 

and is contained within a proper real algebraic (hence measure 0) subset of hA n . 

Thus, we have completed the solution to the problem of how many images of a star can be 
created by a gravitational lens consisting of n point masses. 

1.1. Structure of the Paper. In Section [2l we use basic properties of real algebraic and semi- 
algebraic sets to show that nonsimple polynomials and rational functions lie within proper real 
algebraic subsets of V n , 1Z n , and M. n . The main tool used in the remainder of the paper is the 
extension of the "Argument Principle" to harmonic functions / : C — > C obtained in 0J, which is 
stated precisely in Section |3] We then prove Theorem II . 1 l in Section @J In Section [5J we present 
Rhie's examples lfT5l and then prove Theorem 11.21 We consider the physical case in Section [6l 
presenting a simplified exposition of Petters' lower bound and proving Theorem II .31 

2. Nonsimple Harmonic Functions and Semialgebraic Geometry 

Recall that NV n , N1Z n , and NM. n denote the sets of non-simple polynomials, rational func- 
tions, and "positive massed" rational functions, respectively. 

Proposition 2.1. For any n > 1, 

• NV n is contained in a proper real algebraic subset ofV n , 

• NlZ n is contained in a proper real algebraic subset of1Z n , and 

• NAi n is contained in a proper real algebraic subset of M. n . 

In order to prove Proposition ^. 1 I we will need to take projections of real algebraic sets. However, 
such projections often fail to be real algebraic sets — notable examples include the projection of 
x = y 2 or xy = 1 on to the x-axis. Thus, we will need to work in the realm of semialgebraic 
geometry, see Q] and 0. 
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Definition 2.2. A semialgebraic subset of W l is a finite union of sets given by finitely many poly- 
nomial equations and inequalities with real coefficients. 

We will need the following key properties of semialgebraic sets: 

Tarski-Seidenberg Theorem. Let Abe a semialgebraic subset ofR n+1 and ir : R n+1 — > W 1 , the 

projection on the first n coordinates. Then ft(A) is a semialgebraic subset ofM. n . 

See [fflThm. 2.2.1]. 

Cylindrical Algebraic Decomposition. Any semialgebraic set can be decomposed into finitely 
many sets, each homeomorphic to [0, l] di for some di. 

See flU Thm. 2.3.6]. 

Definition 2.3. The dimension of a semialgebraic set A is the maximum of the dimensions di from 
some cylindrical algebraic decomposition of A. 

Note that the dim(A) is well-defined independent of which cylindrical algebraic decomposition 
is chosen. Recall that for any B C IR n , the Zariski closure of B, denoted Zariski(£>), is the smallest 
real algebraic set containing B. 

Dimension Theorem. For any semialgebraic set A C W" 1 , dim(A) coincides with dim(Zariski(A)). 
See El Section 2.8]. 

Remark. The definition of dimension for a semialgebraic set A is given in a different, but equivalent 
way in 01 Section 2.8]. They define dim(A) = dim (Zariski (A)) and then prove that this definition 
coincides with the maximal dimension of any cell from the Cylindrical Algebraic Decomposition. 

Proof of Proposition ^. 1\ Consider the real algebraic set 



Note that since V is a real algebraic set, it is also a semialgebraic set. 

If we let 7r : V n x C — > V n be the projection onto the first coordinate. Since NV n = 7r(V), the 
Tarski-Seidenberg Theorem gives that NV n is a semialgebraic set. 

We know from [|9] Lemma 5] that SV n is dense in V n , hence NV n cannot contain a set homeo- 
morphic to [0, l] 2n + 2 . Thus, by Definition [23] and the Dimension Theorem, we have that 



In particular, NVn is a subset of Zariski(A r 7 : ' ri ), which is a proper real algebraic subset of V n . 

An identical proof shows that N1Z n and NAi n are semialgebraic subsets of lZ n and Ai n , re- 
spectively. By the lemma from p. 1081 of 0, S1Z n is a dense subset of lZ n , implying that 
dim(Zariski(iV^ n )) < dim(fc n ). 

Moreover, given any r(z) E M. n and any c 6 C, r(z) + c corresponds to a shift made to the 
locations of each of the masses Zi and therefore, r(z) + c G Ai n as well. Thus, the lemma from [|7] 
p. 1081] also shows that SM. n is a dense subset of M. n , implying that dim(Zariski(A r A1 n )) < 
dim(A^ n ). □ 




dim(Zariski(A^P n )) = $xm{NV n ) <2n + l. 
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3. Harmonic Functions and the Argument Principle 

We will need to use an extension of the classical Argument Principle to harmonic functions given 
by flU and |[T8l . For a harmonic function f(x, y) = u(x, y) + iv(x, y) defined on an open simply 
connected set U C C, we can find analytic functions h and g, unique up to additive constants, such 
that f — h + g. Let us consider the power-series expansions of h and g at z : 

oo oo 

h(z) = a + ^2a k (z - z ) k , g(z) = b + ^b k (z - z ) k . 

fc=l fc=l 

Let m > 1 be the first index for which either a m or b m is non-zero. We say that / is sense- 
preserving (s.p) at z if a m ^ and \b m /a m \ < 1, and we say that / is sense-reversing (s.r) at z if 
b m 7^ and \a m /b m \ < 1. Note that if the Jacobian Vf(z ) = \ai\ 2 — \b\\ 2 ^ 0, then this definition 
coincides with the classical one. If h(z ) = 0, we define the order of z as +m if / is s.p at z and 
— m if / is s.r at z . If / is neither s.p or s.r at z , then z is called a singular point and the order is 
undefined. 

We will also need to consider harmonic functions with poles, i.e. functions 

/:C\{zi,. ..,**}-► C, 

which are harmonic and satisfy lini^^ \f(z)\ = oo. The points Zi, . . . , z k are called poles of / 
and we will write f(zj) = oo. 

Take an oriented closed contour T such that f(x,y) {0, oo} for (x,y) E T. Consider a 
"normal" coordinate s : [0, 1] — > T, with s(0) = s(l), and write f(s(t)) = r(t)e l6 ^ in such a way 
that 9 varies continuously over [0, 1]. Then we say that A r arg(/) = #(1) — 0(0) = 2ix ■ u-p, where 
cur is an integer. We call cor the "winding number" of / over T. 

A simple calculation shows that the order of a zero is equal to w 7 , where 7 is a sufficiently small 
geometric circle centered at the zero and positively oriented. The order of a pole is defined in an 
ad-hoc way to be — cu 7 , where 7 is defined similarly. 

The result from |[T8l is: 

Argument Principle. Let f be harmonic, except at a finite number of poles, in a simply connected 
domain D C C. Let C be an oriented closed contour in D not passing through a pole or a zero, 
enclosing a region (taken with orientation) Q C D. Suppose that f has no singular zeros in D, 
and let N be the sum of the orders of the zeros of f in Q. Let M be the sum of the orders of the 
poles of f in Q. Then 

A c arg f(z) = 2ttN - 2ixM. 

4. Polynomial Case 

In this section, we prove Theorem 11.11 An immediate application of the Argument Principle 
gives the following result (briefly mentioned in [|9)). 

Proposition 4.1. Let p be a simple polynomial and let k + and k- denote the number of sense- 
preserving and sense-reversing zeros of f(z) = p(z) — z, respectively. If n — degp > 1, then 
k + — k_ = n. In particular, the total number of zeros satisfies k = k + + k^ = n mod 2. 



Proof Let Sr be a circle sufficiently large and centered at the origin so that Sr contains all zeros 
of / and so that / applied to Sr is dominated by the n-th degree term. Then, by the Argument 
Principle, k + — k_ = u Sr = n. □ 
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Proof of Theorem [7771 We have already shown in Proposition 12. 1 I that NV n is contained in a proper 
real algebraic subset of V n so we will focus our attention on simple polynomials. From Proposi- 
tion |4J3 we have that the number of zeros k = k + + k- > n and that k is congruent to n mod 2. 
From [9], we have that k < 3n — 2. Moreover by the Implicit Function Theorem, the sets SV n (k) 
are open for any k. Thus, it suffices to show that SV n (k) is non-empty for k = n, n + 2, . . . , 3n — 2. 
We will prove this by induction on n > 2. For n = 2, z 2 + 1 is an element of S"P 2 (2) since 
z 2 + 1 = z has solutions z — — | ± ^-i. Also, z 2 is an element of S"P2(4) since z 2 = z has 

solutions z — 0, 1, — | ± ^i. 

Suppose SV n (k) ^ for k = n, n + 2, . . . , 3n — 2. For each such fc, we will show that 
SV n +i( k + 1) ^ 0- We first let 

c(z) 
d(z) 
c(*) 

Since d(z) and c(z) both have degree n — 1, \d(z) \ + |c(^)| < C^l™ -1 for some constant C. Now, 
assume p e SV n (k). Choose i? > sufficiently large such that the following hold true: 

(1) All roots of p(z) — z are contained in the circle Sr of radius R centered at the origin. 

(2) Outside of S R , - C\z\ n ~ l > \z\. 

Denote the k roots of f(z) = p(z) — z as r±, . . . , r^. For each i, let D r . be a sufficiently small disk 
centered at such that the D n are pairwise disjoint and: 

(3) There exists A such that for all z G C\ |J D n , \f(z) \ = \p(z) - z\ > A > 0. 

(4) There exists B such that for all z e{jD n , \V(z)\ = \ \p'(z)\ 2 - 1| > B > 0. 

Let = + e2 n+1 , /(z) = p(z) — z, and let k + and A;_ denote the number of s.p and s.r 
zeros of /, respectively. We will first prove that k_ = k_. Using (1) and (2), we can find e > 
sufficiently small so that 

(5) / remains non-zero inside of and outside of |J D Ti , 

(6) the winding number of / on each of the circles C n = dD r . remains as ±2n, and 

(7) the Jacobians of / and / have the same sign on each D r . . 

Thus, by the Argument Principle, / has exactly one zero within each disc D r . and this zero has the 
same orientation as the original zero of / in D Ti . We conclude that / and / have the same number 
of s.r zeros inside of Sr. 

We now show that / has no s.r zeros outside of Sr. Outside of Sr, f(z) is s.p, except possibly 
in a small region determined by the inequality \j?(z) \ 2 = \(n + l)ez n + na n z n ~ l + c(z)\ 2 < 1. By 
the triangle inequality, all points in the non-s.p region satisfy 

\(n+ l)ez n + na n z n ~ l \ < 1 + \c{z)\. 



p'(z) — na n z n 



{n-l)a n -iz n +... + GH, 

p(z) - a n z n = a n -iz n ~ l + . . . + a , and 

z 



n+1 



\c(z)\). 
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On this region, 



\f(z)\ > \ez n+1 + a n z n \ 



\d(z)\ 



\z\ > 



a n z' 



n 



n 



I (n + l)ez n + na n z 



n-l I 



\d(z)\ 



> 



> 



a n z' 



n + 1 

(l n Z 



n + 1 



n + 1 



(l + \c{z)\)-\d{z)\-\z\> 



a n z' 



n 



\c{z)\-\d{z)\-\z\ 



\z\. 



By the assumption on R, \f(z)\ > on the non-s.p region outside of Sr. Since / and / have 
the same number of s.r zeros inside of Sr and neither of them has any zeros outside of Sr, we 
conclude that A;_ = . 

By Proposition 14. 1[ k + — = n + 1, so k + = k + + 1, which gives that k = k + 1. Thus we 
have shown that SV n +i(k + 1) ^ for k + 1 = n + 1, n + 3, . . . , 3n — 1. It remains to show 
SV n +i(3(n + 1) — 2) 0; however, this follows from Geyer's examples [|6]. □ 



5. Rational Case 

In this section, we will prove Theorem 11.21 using methods similar to those in the polynomial 
case. Throughout the proof, we will only consider rational functions r(z) = |?4 G 1Z n for which 
deg q > deg p. In terms of the description of lZ n given in the introduction, this amounts to throwing 
out the algebraic hyperplane b n = 0; see Equation [L3l 

Using the Argument Principle, we can prove the following lemma, which has previously ap- 
peared in [7 , Cor. 2]. 

Lemma 5.1. Let r{z) = be a simple rational function with deg q > degp and let k + and 

denote the number of sense-preserving and sense-reversing roots of f(z) = r(z) — z, respectively. 
Then k + — = n — 1. 

Proof of Lemma \5J\ Let f(z) = r(z) — z. Consider a circle Sr of large radius R centered at 
the origin such that the Sr contains all the zeros and poles of /. Since degp < degg, r(z) is at 
most 0(1) for z large. Then / is sense-reversing on Sr with an argument change of — 2n. By the 
Argument Principle, — 2n = 2n ■ (N — M), where N = k + — k- and M is the number of poles of / 
counted with orientation. But notice that at each of the n = deg q poles of r, / is sense-preserving, 
so — 2n = 2n ■ (k + — k_ — n), which gives k + — k_ — n — 1, as desired. □ 

The examples given in Proposition [5T2] below were previously presented by Rhie in a preprint 
lfT31 . Since they were never published and since we will later use details from the construction, we 
will reproduce them here. 

Proposition 5.2 (Rhie lfT5lO . For n > 2, SlZ n {hn — 5) is a non-empty, open subset oflZ n . 

Proof. We will use the Lens Equation (11.21) to show that the upper bound of 5n — 5 solutions is 
attainable for n > 2 by a specific choice of masses o~j and locations of masses Zj. Throughout the 
proof, the reader may find it helpful to look at Figure |2] where the construction is illustrated with 
four and five masses. 

Consider n equal masses (with one mass on the positive real axis) equally spaced around a circle 
of radius a = (n— 1)~« 2 ( m f act > an Y small a will work) centered at 0. With Zj = ae 2 ^- 7 " 1 ), 



COUNTING ZEROS OF HARMONIC RATIONAL FUNCTIONS 



9 




Figure 2. The diagram on the left shows the real and imaginary parts of the lens 
equation for three point masses around a circle and an epsilon mass at the center. 
The diagram on the right shows the configuration for four point masses and an 
epsilon mass at the center. The bold M's represent mass positions. The black dots 
are sense-preserving image positions and the gray dots are sense-reversing image 
positions. 



(11.21) becomes: 

z n ~ 1 

= r(z) — z; r(z) = 



In the case that n is odd, we only need to consider the lens equation on the real line since the 
mass configuration is symmetric by a rotation of — . There is always one solution at z = 0, and 
the other solutions satisfy 

z n - z n ~ 2 - a n = 0. (5.1) 

The equation z n —z n ~ 2 has three real solutions, and it is easy to see that shifting this graph vertically 
by a sufficiently small amount (a n in our case) will still give three real solutions. Hence, we have 
a total of 3n + 1 images. (This was previously observed in [[TOl .) 

In the case that n is even, we need to consider the real line and the line obtained by a rotation of 
ir/n. Solving (15.11) gives two real solutions by an analysis similar to that in the odd case. On the 
rotated line, we set z = te t7T ^ n with t ^ and the lens equation, 

z n-x 



simplifies to 

t n _ t n-2 + fl n = 0. (5.2) 
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A similar analysis of this equation gives four real solutions. Therefore, there are a total of six 
non-zero solutions on the two neighboring lines, and the total number of images in the even case 
is6-ra/2+l = 3n+l. 

We now drop a small mass e > at the center of the circle. This changes the number of masses 
to n + 1, so we now require 5n distinct image positions. We will first check that if e is sufficiently 
small, then all 3n solutions away from the origin will persist and any new solutions will occur 
within distance 2y/e from the origin. 

Let f(z) = r(z) — z. Denote the 3n + 1 roots of / as = r , . . . , r 3n . For each i, let D n be a 
sufficiently small disk centered at r.j such that the D r . are pairwise disjoint and 

(1) There exists A such that for all z e C \ |J D n , \f(z)\ = \r(z) - z\ > A > 0. 

(2) There exists B such that for all z £ \JD n , \T>(z)\ = \ \r'{z) 2 \ - 1| > B > 0. 

Suppose that e is sufficiently small so that A > sfe and so that the disk D' of radius 2^/e centered 
at satisfies D' rQ C D ro . By (2), there are no critical points in the annulus D ro \ D' rQ , so the 
minimum of |/| (over this annulus) occurs on its boundary. On dD ro , \f(z)\ > B, and on dD' 

\f{z)\ = \r{z)-z\>\z\-\r{z)\>y/l 
since r(0) = r'(0) = 0. Hence, we have: 

(10 For all z G C \ D' ro U D n U ■ • • U D r3n , \f(z)\ > sfe. 

Adding the mass e results in f(z) = r(z) — z = r(z) + - — z. By (10 and (2), if e is sufficiently 

small, then in the region \z\> 2y/e, each of the zeros of / is obtained by a continuous motion of a 
zero of /, and each of them will remain simple with the same orientation. 
We now focus on\z\ < 2y/e. The new lens equation on the real line is 

z n+2 _ z n^ z 2 a n + £fl n = Q (53) 

At this scale, equation (15.31) is an arbitrarily small perturbation of the equation — z 2 a n + ea n = 0, 
which gives two simple solutions at ±y/e. If n is odd, then symmetry under rotation by 2n /n leads 
to a total of 2n new simple images in the 2y/e neighborhood of the origin. Since we originally had 
3n + 1 simple zeros, we now have exactly exactly 3n + 1 + 2n — 1 = 5n simple zeros, as desired. 

If n is even, symmetry under rotation by 2n /n only leads to n new images in the 2^/e neigh- 
borhood of the origin. However, on the line z = te ln ^ n the lens equation is the same as Equation 
(15.31) . except that the first two terms have the opposite sign. Two new simple zeros are produced 
on this line at approximately t = ±y/e. Under symmetry by rotation of 2n/n, this produces the n 
additional new images that we needed. 

In the case of n = 2 and n = 3, this construction fails due to the low number of masses. 
However, a quick check shows that the original construction of n equally spaced masses around 
the circle of radius a gives 5 and 10 simple zeros in these two cases. For example, in the case of 
n = 2, we have r(z) = z il x i 2 -> an d solving the equation z = r(z) will give five solutions. The 
reason we do not reach the 3n + 1 bound in the case n = 2 is because on the imaginary axis, there 
are only two solutions (since the degree is not high enough). In the case of n = 3, 3n + 1 = 5n — 5, 
so we see that the original construction works as well. □ 

Proposition 5.3. For n > 2, S1Z n (5n — 7) is a non-empty, open subset oflZ n . 

Proof of Proposition ^. 3\ Placing mass 1/2 at z = 1 and z — — 1 produces three images. Placing 
mass 1 at z — \ and z = — \ ± produces eight images. 

For four or more masses, the construction is a perturbation of Rhie's examples, which consisted 
of n equal masses equally spaced around the circle of radius a and an e mass at the origin. We 
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claim that a small perturbation for the location of the e mass from the origin will result in 5n — 2 
simple solutions. More exactly, for appropriate choices of A and 0, we claim that the equation 

n-l e 

+ 7 = z, (5.4) 



where b = Ae( n ~ 1 " 2 ■ e 1 ^ and c = a n is the constant from Proposition 15 .2[ has exactly hn — 2 
simple solutions so long as e > is sufficiently small. 

As in the Proof of Proposition |5?2l a perturbation by -§r of f(z) = r(z) — z will not affect the 
roots outside of the disk of radius 2 y/e centered at the origin so long as \b\ < y/e. Hence, we need 
only consider Equation (15.41 ) inside the disk of radius 2 y/e. 

Equation (15.41) . after simplification, becomes 

z n , bz n ~ l 
h e = zz - bz H . (5.5) 



Replacing z n — c by — c in Equation (15.51 ) gives 

z n bz 71 ^ 1 
e = zz — bz , (5.6) 

c c 

which is an approximation of Equation (15.51 ) up to an error term of order e n . Let z = re %e . The real 
part of Equation (15.61 ) gives r = y/e + 0(e). We now compare the imaginary parts of both sides, 
which will give us the values of 9. 

If n is even, let = 7r/2 + 7r/2n. The lowest order terms for the imaginary part of Equation (15.61) 
are of order e n l 2 . They are 

Imf-^ = \m(-bz). 

Using that r = y/e + 0(e) and |6| = Ae~ n ~ 1 " 2 , the lowest order terms give 

sinn^ = C sin(0 — 9), 

where C = A ■ c. We notice that when 9 = ir/2n, both sides attain their maximum value. The left 
hand side attains its minimum value when 9 = g*±ij)ZE f or integer k. The right hand side attains 
its minimum when 9 = (~ 4fcn ~^ n+1 ) 7r f or m t e g er ^ However, 4k + 3 ^ —Akn — 2n + 1 (mod 4) 
since n is even. Therefore, sinn9 and sin(0 — 9) are tangent at only one point, 9 = n/2n. When 
C = 0, there are 2n solutions - a simple calculation shows that for C < 1, sinn^ and sin(0 — #) 
cannot be tangent. Hence, the first time a tangency occurs is when C = 1, and from above, this 
tangency is located at the single point 9 = n /2n. Increasing C by a sufficiently small amount will 
pull apart this tangency, but the other simple solutions to sinn# = sin(0 — 9) will persist, so we 
have lost exactly two solutions, as desired. 

When n is odd, for any choice of 0, the above approach produces two tangencies and thus it 
cannot be used without considering higher order terms. Let 4> = n/2. As before, r = y/e + 0(e). 
The imaginary part of Equation (15.61 ) yields 

sinn# = ^Lcos# + — cos((n - 1)9). (5.7) 

Up to terms of order e^ n_1 ^ 2 , Equation (15.71 ) becomes 

sinn9 = Ccos9 + Ae {n ~ 2)/2 cos((ra - 1)0), 
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where C = A ■ c. We first look at the lowest order terms, which gives 

sin ri0 = Ccos0. (5.8) 

When the mass e is located exactly at the center, C — 0, and there are 2n solutions. First, notice 
that if 9 is a solution to (15.81) . so is ir + 9. Increasing C will not change the number of solutions 
to (15.81) until C cos 9 is tangent to sin n9 at some point. For this to occur, we require a solution to 
(15.81) to also be a solution to 

n cos n9 = — Csin#. (5.9) 
Squaring ra- (!5.8l) and (15.91 ) and summing gives 

n 2 = C 2 + C 2 (n 2 -l)-cos 2 9, 

which gives two solutions for cos 9 and four for 9 - however, one can readily verify that two of 
these are extraneous, so there are only two values of 9 for which (15.81 ) and (15.91 ) are satisfied, and 
moreover they differ by n. Reconsidering the small term Ae^ n ~ 2 ^ 2 cos((n — 1)9), we find that 
cos((n — 1)9) = cos((n — l)(ir + 9)). Since one of the tangencies was a minimum for both sides of 
Equation (15.81) and the other was a maximum, this perturbation will pull apart one of the tangencies 
and change the other to be two simple solutions. □ 

With these facts, we can prove Theorem II .21 

Proof of Theorem \1.2\ We have already shown in Proposition 12. 1 I that NlZ n is contained in a proper 
real algebraic subset of lZ n so we will focus our attention on simple rational functions. 

We will also throw out the proper algebraic set given by b n = in the parameterization of 
rational functions (Equation (II .3b ). This allows us to restrict our attention to the case that deg q > 
degp. Under these assumptions, we have that the number of zeros, k, satisfies n — 1 < k, from 
the Argument Principle, and k < 5n — 5, from [|9]. Furthermore, from Lemma ISTTl we have 
that k = n — 1 mod 2. Thus, it suffices to show that for each n > 2, SlZ n {k) is non-empty for 
k = n — l,n + l,7i + 3, ...,5n — 5. 

We will proceed with an induction on n, similar to the proof of Theorem 11.11 First consider 
the case where n = 2. The equation r(z) = -% = z has exactly one root for any c so STZ 2 (1) is 
non-empty. Furthermore, by Proposition 15 .2 l and Proposition l5.3l we have that S1Z2(3) and STZ2(5) 
are non-empty, respectively. 

Now, suppose for some n > 2 that S1Z n (k) ^ for k = n — 1, n + 1, n + 3, . . . , 5n — 5. 
We will first show that SlZ n+ i(k + 1) ^ for the same values of k. By assumption, there exists 
r(z) E S1Z n (k). Denote the k roots of f(z) = r(z) — z as r±, . . . , r k and the n poles of / as 
zi, Z2,..., z n . Consider the rational function r = r(z) + -+=, where z ^ rj and z ^ z t for any i. 
For each i, let D n be a sufficiently small disk centered at and let D Zi be a sufficiently small disk 
centered at Zi such that the D n and D Zi are pairwise disjoint, z |J D n |J D Zi , and 

(1) There exists A such that for all z G C\[jD ri , \f(z)\ = \z - r(z)\ > A > 0. 

(2) There exists B such that for all z e[jD n , \V(z)\ = \ \r'(z) 2 \ - 1| > B > 0. 

(3) There exists C such that for all z £ C\ \jD Zi , \r'(z)\ < C. 

Clearly, r has degree n+ 1, so it is sufficient to show that / = r(z) — z has k + 1 zeros. For e > 
sufficiently small, the winding number around each of the circles C Ti = dD r . remains as ±27r and 
the Jacobian does not change sign in D r .. Therefore, each root of / moves continuously to some 
new simple root of the same "orientation" as r^. Moreover, Fj remains in the disk D n , so there 
is still exactly one root in each disk. 
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We claim now that any new root of / must be sense-preserving. If f(z) = 0, then 

A 2 



z — z 



At any such point, 



1^0)1 



r'{z) 



{z-z)' 



r(z)\ > A 



> 



(z-z)' 



> 



(z-I)< 



A 2 

- \r'(z)\ > — -C. 



As long as e is sufficiently small, \r'(z)\ > 1, so any new root of / is sense-preserving. By Lemma 
15.11 k + — k_ = n. As no new s.r roots are created, stays the same, so k + must have increased 
by exactly 1, as desired. 

We have shown that S n+ i(k) is non-empty for A; = n, n + 2, . . . , 5n — 4 = 5(n+l) — 9. However, 
Propositions I5.2l and l5.3l give that S n+ x(5n — 2) and S n+ i(bn) are also non-empty, thus our proof 
is complete. □ 

6. Physical Case 



The Lens Equation (11.21) can be rewritten as 







Oi(z - Zi) 



~ 12 



(6.1) 



where is the source position, Zj is the position of the ith point mass, and z is the image position. 
The time delay function T is defined to be 

T(z) = — y^Ojln|z - Zi\. 



Note that zq = x + iy is a solution to d z T = (T x — iT y )/2 = (a critical point of T) if and only 
if z is a solution to Equation 16.11 Moreover, z is a non-degenerate local minimum (respectively 
maximum) of T iff z is a s.r (respectively s.p) simple zero of (16. It . 

Proposition 6.1 (Petters [[ill ). If all of the masses cr; are positive and if all of the zeros of the lens 
equation (16.11) are simple, then there are at least n + 1 of them. 

Proof. Let D R be the closed disc centered at the origin of radius R and for each i, let D ij6 be the 
open disc of radius e centered at Zj. Let G = D R \ UD i e . 

If | z | —7- oo or if z — > Zj, then T(z) — > oo. Hence, we can choose i? > sufficiently large and 
e > sufficiently small so that the minimum of T is not attained on the boundary of G. Therefore, 
the minimum of T occurs on the interior of G. This minimum corresponds to a s.r simple zero of 
( 16.11) . We can show that n — 1 = k + — as in Lemma [5TT1 where k + (&_) is the number of s.p (s.r) 
zeros of (I6.ll ). Since we have at least one s.r zero, k_ > 1, and k + + > n — 1 + 2fc_ > n + 1, 
as desired. □ 



Proof of Theorem [L3\ We use an inductive proof similar to that of Theorem II .21 

Consider the case n = 2. Placing mass 1 at any point other than the origin will result in two 
images and the inductive step from the proof of Theorem II .2 I shows that placing a sufficiently small 
e mass anywhere else will generate exactly one more image. Also, Rhie lfl"5l gives an example for 
a configuration that yields 5 roots. Thus, S'A^2(3) and SA^^) are non-empty open sets. 

From here, exactly the same inductive step can be applied as in the proof of Theorem II .2 I to show 
that if SMniji + 1), . . . , SM n (5n - 5) are non-empty, then SM n+1 (n + 3), . . . , SM n+1 (5n - 4) 
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are non-empty since addition of the term e/(z — z) corresponds to adding a small mass at z. 
Propositions 15 .21 and 15 .3 1 give that SM n+ i(hn — 2) and SM n+ i(5n) are non-empty. □ 

Remark. Each of the SM. n {k) contains a positive massed rational function with masses arbitrarily 
close to the origin. Thus, these functions are consistent with the physical assumption that all the 
masses are close together and at small angular positions with respect to the light source. In the 
proof of Theorem II .31 the base case of the induction can be done with two masses arbitrarily close 
to the origin. The inductive step from the proof of Theorem II . 21 allowed us to add a small mass at 
an arbitrary point and, moreover, the proofs of Propositions 15.21 and 15.31 can be done with masses 
arbitrarily close to the origin. 
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